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Abstract 

In this work we prove that the unique 1-convex solution of the Monge- 
Kantorovitch measure transportation problem between the Wiener mea- 
sure and a target measure which has an //-log-concave density, in the sense 
of [Sj, w.r.to the Wiener measure is also the strong solution of the Monge- 
Ampere equation in the frame of infinite dimensional Frechet spaces. We 
further enhance the polar factorization results of the mappings which trans- 
form a spread measure to another one in terms of the measure transporta- 
tion of Monge-Kantorovitch and clarify the relation between this concept 
and the Ito-solutions of the Monge- Ampere equation. 

1 Introduction 

In 1781, G. Monge has launched his famous problem ^7], which can be expressed 
in terms of the modern mathematics as follows: given two probability measures 
p and v on IR n , find the map T : IR n — > IR n such that Tp = v 1 and T is also the 
solution of the minimization problem 

inf i / c(x, U{x))p(dx) \ , (1.1) 



u 



IR" 



where the infimum is taken between all the maps U : H n — > lR n such that Up = v 
and where c : H n x IR n — > IR + is a positive, measurable function, called usually 
the cost function. In the original problem of Monge, the cost function c(x, y) was 
| a; — y\ and the dimension n was three. Later other costs have been considered, 
between them, the most popular one which is also abundantly studied, is the 
case where c(x,y) — \x — y\ 2 . After several tentatives (cf., [UI2]), m the 1940's 
this highly nonlinear problem of Monge has been reduced to a linear problem by 

1 Tp means the image of the measure p under the map T 
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Kantorovitch, cf.^3], in the following way: let E(p, v) be the set of probability 
measures on H n x ]R n , whose first marginals are p and the second marginals are 
v. Find the element (s) of E(p, v) which are the solutions of the minimization 
problem: 

inf ( / c(x,y)d(3(x,y)\ . (1.2) 

It is obviuous that £(p, v) is a convex, compact set under the weak*-topology 
of measures, hence, in case, the cost function c has some regularity properties, 
like being lower semi-continuous, this problem would have solutions. If any one 
of them is supported by the graph of a map T : ]R n — > H ra , then obviously, T 
will be also a solution of the original problem of Monge [T7T1 Since that time, the 
problem (jl.2j) is called the Monge-Kantorovitch problem (MKP) . The program of 
Kantorovitch has been followed by several people and a major contribution has 
been done by Sudakov [22,. In the early 90's there has been another impetus to 
this problem, cf., [HUH], where it has been discovered the important role played 
by the convex functions in the construction of the solutions of the MKP and of 
the problem of Monge (cf., [1^1 HE1 HE] ) . We refer the reader to jH] and to [2H] for 
recent surveys. 

In [TT], we have solved the MKP and the problem of Monge in the infinite 
dimensional case, where the measures are concentrated in a Frechet space W into 
which a Hilbert space H is injected densely and continuously. We call H the 
Cameron-Martin space in reference to the Gaussian case. The cost function is 
defined on W x W as 

c(x,y) = \x — y\ 2 H if x — y E H 
= oo if x — y j£ H , 

where | • \h denotes the Euclidean norm of H. Because of this choice, in comparison 
to the finite dimensional space, the situation becomes quite singular, since, in 
general, the Cameron-Martin space if is a negligeable set (i.e., of null measure) 
with respect to almost all reasonable measures for which one can expect to have 
solutions of the problems of Monge and of MKP. On the other hand, due to 
the potential applications to several problems of analysis and physics, this cost 
function is particularly important. For example, it is particularly well- adapted 
to the study of the absolute continuity of the image of the Wiener measure under 
the perturbations of identity, which is a subject under investigation since the 
early works of N. Wiener, R.H. Cameron and W.T. Martin and of several other 
mathematicians and engineers who have made worthy contributions (cf. the list 
of references of |26j). 

This paper is devoted to the further developments of the subject. At first 
we give a generalization of the polar factorization of vector fields which map a 
probability measure on W to another one such that one of them is spread (cf. the 
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preliminaries) and the two measures are at finite Wasserstein distance from each 
other (without any absolute continuity hypothesis). As an example we treat in 
detail the case of the infinite dimensional Gaussian measures. 

The proof of the fact that the transport map, when the target measure has 
an if -log-concave density, satisfies the functional analytic (or strong) Monge- 
Ampere equation is probably the most important contribution of this paper. 
In [TT], we have studied the Monge- Ampere equation for the upper and lower 
bounded densities with respect to the (infinite dimensional) Wiener measure. 
The main difficulty in this infinite dimensional case stems from the lack of reg- 
ularity of the transport potentials, in fact we only know that these functions 
are in the Gaussian Sobolev space ID 2i i, i.e., they have only first order Sobolev 
derivatives. However, to write the Gaussian Jacobian, we need them to have 
second order Sobolev derivatives taking values in the space of Hilbert-Schmidt 
operators on the Cameron-Martin space H. This difficulty is worse than those we 
encounter in the finite dimensional case, since in the latter the Hilbert-Schmidt 
property holds automatically. Moreover, in the finite dimensional situation the 
lack of second order derivatives is solved with the help of the Alexandroff deriva- 
tives of the convex functions. In the infinite dimensional case the situation is 
worse: the transport potentials are not in general convex, nor H-convex (which 
is a more reasonable requirement than being convex, cf. 9\), but only 1-convex 
in the Cameron-Martin space direction. Hence their second order derivatives in 
the sense of distributions are not in general measures; even if this happens in 
some exceptional situations, their absolutely continuous parts do not take values 
in the space of Hilbert-Schmidt operators, a condition which is indispensable to 
write down the Jacobian of the transport map. Hence it is impossible in gen- 
eral to construct the strong solutions of the Monge- Ampere equation. In Section 
combining the finite dimensional results of Caffarelli [7j with Wiener space 
analysis, we solve completely this problem when the target measure is if-log- 
concave. More precisely, we show that the transport potential has a second order 
derivative as an Hilbert-Schmidt operator valued map, hence we can write the 
corresponding Jacobian which includes the modified Carleman-Fredholm deter- 
minant, cf. ,26^ and finally we prove that the transport potential is the unique 
1-convex strong solution of the Monge- Ampere equation. In Section El we show 
that all these difficulties disappear if we use the natural Ito Calculus and we can 
calculate the Ito Jacobian (cf. Theorem I6.1J1 using the natural Brownian mo- 
tion which is associated to the solution of the Monge problem. In fact, with Ito 
parametrization, the complications are absorbed by the nitrations of forward and 
backward transport processes (i.e., maps). We give also the delicate relations be- 
tween the polar factorization of the absolutely continuous transformations of the 
Wiener measure and the Brownian motions which appear in the semimartingale 
decomposition of the transport process with respect to its natural filtration. 
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2 Preliminaries and notations 

Let W be a separable Frechet space equipped with a Gaussian measure \i of zero 
mean whose support is the whole space 2 . The corresponding Cameron-Martin 
space is denoted by H. Recall that the injection H W is compact and its 
adjoint is the natural injection W* c — > H* C L 2 (fi). The triple (W,fi,H) is 
called an abstract Wiener space. Recall that W — H if and only if W is finite 
dimensional. A subspace F of H is called regular if the corresponding orthogonal 
projection has a continuous extension to W, denoted again by the same letter. 
It is well-known that there exists an increasing sequence of regular subspaces 
(F n ,n > 1), called total, such that U n F n is dense in H and in W. Let V n be 
the cr-algebra generated by 7i> n , then for any / G L p (fi), the martingale sequence 
(E[f\V n ],n > 1) converges to / (strongly if p < oo) in L p (fi). Observe that the 
function /„ = i£[f|Vy can be identified with a function on the finite dimensional 
abstract Wiener space (F n , /i n , F n ), where fi n = Tr n fi. 

Since the translations of \i with the elements of H induce measures equivalent 
to fi, the Gateaux derivative in H direction of the random variables is a closable 
operator on L p (/i)-spaces and this closure will be denoted by V cf., for exam- 
ple [23] • The corresponding Sobolev spaces (the equivalence classes) of the real 
random variables will be denoted as HDj,*,, where k G IN is the order of differen- 
tiability and p > 1 is the order of integrability. If the random variables are with 
values in some separable Hilbert space, say then we shall define similarly the 
corresponding Sobolev spaces and they are denoted as E) p fc ( ( l ) ), p > 1, k 6 IN. 
Since V : TD P ± — » ^D P ,k-i(H) is a continuous and linear operator its adjoint is 
a well-defined operator which we represent by 5. In the case of classical Wiener 
space, i.e., when W = C(lR + ,IR d ), then 5 coincides with the It 6 integral of the 
Lebesgue density of the adapted elements of JD Pt k(H) (cf.[25j). 

For any t > and measurable / : W — > 1R+, we note by 



it is well-known that (Pt,t G IR+) is a hypercontractive semigroup on L p (fj,),p > 
1, which is called the Ornstein-Uhlenbeck semigroup (cf.[25|). Its infinitesimal 
generator is denoted by — C and we call C the Ornstein-Uhlenbeck operator (some- 
times called the number operator by the physicists). Due to the Meyer inequalities 
(cf., for instance |25j). the norms defined by 



are equivalent to the norms defined by the iterates of the Sobolev derivative 
V. This observation permits us to identify the duals of the space ID P) fe($);p > 

2 The reader may assume that W = C(M+, H d ), d > 1 or W = IR^. 




f\\p,k 



||(/ + £) fe/2 



(2.3) 
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1, k £ IN by Dg 5 _fc($'), with q^ 1 = 1 — p^ 1 , where the latter space is defined by 
replacing k in ()2.3|) by —k, this gives us the distribution spaces on the Wiener 
space W (in fact we can take as k any real number). An easy calculation shows 
that, formally, 5 o V = £, and this permits us to extend the divergence and the 
derivative operators to the distributions as linear, continuous operators. In fact 
5 : TD q , k (H ® $) -> H\ fe _i($) and V : B, iA .($) -> E) g)fc _i(if ® $) continuously, 
for any q > 1 and £ H, where if Cg> $ denotes the completed Hilbert-Schmidt 
tensor product (cf., for instance j2H])- The following assertion is useful: assume 
that (Z n ,n > 1) C ID' converges to Z in ID', assume further that each each Z n 
is a probability measure on W, then Z is also a probability and (Z n ,n > 1) 
converges to Z in the weak topology of measures. In particular, a lower bounded 
distribution (in the sense that there exists a constant c £ 1R such that Z + c is a 
positive distribution) is a (Radon) measure on W, cf. 

A measurable function / : W — > IR U {oo} is called if-convex (cf.[9j) if 

/i -> f(x + h) 

is convex //-almost surely, i.e., if for any h,k £ if, s,t £ [0, 1], s + £ = 1, we have 

f(x + sh + tk) < sf(x + h)+ tf(x + k) , 

almost surely, where the negligeable set on which this inequality fails may depend 
on the choice of s, h and of k. We can rephrase this property by saying that 
h — > (x — > f(x + h)) is an L°(/i)-valued convex function on H. f is called 
1-convex if the map 

h^(x^f(x + h)+±\h\%\ 

is convex on the Cameron-Martin space H with values in L°(n). Note that all 
these notions are compatible with the /i-equivalence classes of random variables 
thanks to the Cameron-Martin theorem. It is proven in |Uj that this definition 
is equivalent the following condition: Let (7r n ,n > 1) be a sequence of regular, 
finite dimensional, orthogonal projections of H, increasing to the identity map 
In- Denote also by 7r n its continuous extension to W and define 7r^ = Iw — ftn- 
For x £ W, let x n = n n x and x^ = n^x. Then / is 1-convex if and only if 

%n * g I ^ n I H f (*^n ) 

is 7r^-/i-almost surely convex. We define similarly the notion of if -concave and 
if -log- concave functions. In particular, one can prove that, for any if -log-concave 
function / on W, P t f and £7[/|V^] are again if- log-concave 0. 

3 Monge-Kantorovitch problem 

Let us recall the definition of the Monge-Kantorovitch problem in our case: 
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Definition 3.1 Let p and v be two probability measures on W , let also £(p, v) be 
the convex subset of the probability measures on the product space WxW whose 
first marginal is p and the second one is v. The Monge-Kantorovitch problem 
for the couple (p,v) consists of finding a measure 7 G £(p, v) which realizes the 
following infimum: 



The function c(x,y) = \x — y\ 2 H is called the cost function. 

Remark 3.1 Note that the cost function is not continuous with respect to the 
product topology of W x W and it takes the value 00 very often for the most 
notable measures, e.g., when p and v are absolutely continuous with respect to 
the Wiener measure p. 

The proof of the next theorem, for which we refer the reader to ^T] , can be done 
by choosing a proper disintegration of any optimal measure in such a way that 
the elements of this disintegration are the solutions of finite dimensional Monge- 
Kantorovitch problems. The latter is proven with the help of the measurable 
section-selection theorem. 

Theorem 3.1 (General case) Suppose that p and v are two probability mea- 
sures on W such that 



Let (7f n ,n > 1) be a total increasing sequence of regular projections (of H , con- 
verging to the identity map of H). Suppose that, for any n > 1, the regular 
conditional probabilities p(- \tt^ = x^) vanish n^p-almost surely on the subsets of 
(7r^")~ 1 (W / ) with Hausdorff dimension n—1. Then there exists a unique solution of 
the Monge-Kantorovitch problem, denoted byy G X(p, v) and'j is supported by the 
graph of a Borel map T which is the solution of the Monge problem. T :W -^-W 
is of the form T = Iw + £ , where £ G H almost surely. Besides we have 



and for n^p-almost almost all x^, the map u — > u + £(u + x^) is cyclically 
monotone on (7r^) _1 {^} ; in the sense that 




inf / 




d H (p, v) < 00 . 




N 




i=l 
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Ti^p-almost surely, for any cyclic sequence {u%, . . . , un,un+i = wi} from ir n (W). 
Finally, if, for any n > 1, TT^u-almost surely, v[- |-7r„ = y L ) also vanishes on the 
n — l-Hausdorff dimensional subsets of {tt^)~ 1 {W) , thenT is invertible, i.e, there 
exists S : W — > W of the form S = Iw + V such that rj e H satisfies a similar 
cyclic monotononicity property as £ and that 



Remark 3.2 In particular, for all the measures p which are absolutely continu- 
ous with respect to the Wiener measure p, the second hypothesis is satisfied, i.e., 
the measure p(- \tc^ = x^) vanishes on the sets of Hausdorff dimension n — 1. 

Any probability measure satisfying the hypothesis of Theorem 13.11 is called a 
spread measure. Namely, 

Definition 3.2 A probability measure m on (W, B(W)) is called a spread mea- 
sure if there exists a sequence of finite dimensional regular projections (ir n , n > 1) 
converging to Ih such that the regular conditional probabilities m( ■ = x^) 
which are concentrated in the n-dimensional spaces n n (W) + x^ vanish on the 
sets of Hausdorff dimension n — 1 for 7r„ (m)- almost all x^ and for any n > 1. 

Definition 3.3 Let m be a probability on (W, B(W)) and g be a real-valued ran- 
dom variable. For any a G [0, 1], we say that g is partially a-convex m-a.s., if 
for any regular finite dimensional subspace F of W , the partial map 



is convex m(-|7r^ = yp)-almost surely, where up is the projection corresponding 
to F, nj; = I w - ti f , y F = ir F (y) and yp- = y - y F . 

The case where one of the measures is the Wiener measure and the other is 
absolutely continuous with respect to p is the most important one for the ap- 
plications. Consequently we give the related results separately in the following 
theorem where the tools of the Malliavin calculus give more information about 
the maps £ and rj of Theorem 13.11 



1 



j{(x,y)eWx W:ToS(y)=y} 
7 {(x, y) e W x W : S o T(x) = x} . 



In particular we have 




yF -> ^Mh + g{yF + y F ) 
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Theorem 3.2 (Gaussian case) Let v be the measure dv = Ldfi, where L is 
a positive random variable, with E[L] = 1. Assume that <i#(/i, v) < oo (for 
instance L G ILloglLj. Then there exists a 1- convex function ip G E>2,i and a 
partially 1- convex function if) G L 2 [y), both are unique upto a constant and called, 
respectively, forward and backward Monge-Kantorovitch potentials, such that 

cp(x) +ip{y) + -y\ 2 H >0 

for all (x, y) G W x W and that 

^(x)+ip(y) + -\x-y\ 2 H = 

^-almost everywhere. The map T = I\y + Vip is the unique solution of the 
original problem of Monge. Moreover, its graph supports the unique solution of 
the Monge-Kantorovitch problem 7. Consequently 

(I w x T)fi = 7 

In particular T maps \x to v and T is almost surely invertible, i.e., there exists 
some T~ l = Iw + rj such that T~ x v = [i, r\ G L 2 (u) and that 

1 = \i {x : T _1 o T(x) = x} 

= u{yeW : ToT-\y) = y) . 

Remark 3.3 Assume that the operator V is closable with respect to v , then we 
have rj = Vip. In particular, if v and \i are equivalent, then we have 

T'^Iw + Vif), 

where ip is a 1-convex function. 

Remark 3.4 Let (e n ,n G IN) be a complete, orthonormal in H, denote by V n 
the sigma algebra generated by {Sex, . . . , d~e n } and let L n = E[L\V n }. If <p n G 1)2,1 
is the function constructed in Theorem W.ty corresponding to L n , then, using the 
inequality (cf., TTTp 

d 2 H {w)<2E[L\ogL], 
we can prove that the sequence ((p n ,n G IN) converges to ip in D 2i i. 
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4 Polar factorization of mappings between spread 
measures 

In pi] we have proved the polar factorization of the mappings U : W — > W such 
that the Wasserstein distance between U\i and the Wiener measure /i, denoted 
by dn(p,Up), is finite. We have also studied the particular case where U is a 
perturbations of identity, i.e., it is the form Iw + u, where u maps W to the 
Cameron-Martin space H . In this section we shall generalize this results in the 
frame of spread measures. 

Theorem 4.1 Assume that p and v are spread measures with dn{p, v) < oo and 
that Up = v, for some measurable map U : W — > W. Let T be the optimal trans- 
port map sending p to v, whose existence and uniqueness is proven in Theorem 
\3.1\ Then R = T~ 1 oU is a p-rotation (i.e., Rp = p) and U = ToR ) morover, if 
U is a perturbation of identity, then R is also a perturbation of identity. In both 
cases, R is the p-almost everywhere unique minimal p-rotation in the sense that 



[ \U(x) - R(x)\ 2 H dp(x) = inf f \U(x) - R'(x)\ 2 H dp(x) 



(4.4) 



where TZ denotes the set of p-rotations. 

Proof: Let T be the optimal transportation of p to v whose existence and unique- 
ness follows from Theorem 13. II The unique solution 7 of the Monge-Kantorovitch 
problem for E(p, v) can be written as 7 = (J x T)p. Since v is spread, T is in- 
vertible on the support of v and we have also 7 = (T _1 x T)v. In particular 
Rp = T _1 o Up — T~ l v = p, hence R is a rotation. Let R' be another rotation in 
1Z, define 7' = (R' x U)p, then 7' e S(p, v) and the optimality of 7 implies that 
^(7) < besides we have 

/ \U(x)-R(x)\ 2 H dp(x) = [ \U(x)-T- 1 oU(x)\ 2 H dp(x) 
Jw Jw 

= [ \x-T-\x)\ 2 H dv(x) 
Jw 

= I \T(x) — x\ 2 H dp{x) 
Jw 

= J(l)- 

On the other hand 

J( 7 ')= / \U(x)-R'(x)\ 2 H dp(x), 
Jw 

hence the relation f)4.4j) follows. Assume now that for the second rotation R' G 1Z 
we have the equality 

\U{x) - R{x)\ 2 H dp{x) = [ \U(x) - R'(x)\ 2 H dp(x). 
w Jw 
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Then we have J(7) = J(j'), where 7' is defined above. By the uniqueness of the 
solution of Monge-Kantorovitch problem due to Theorem 13.11 we should have 
7 = 7'. Hence (R x U)p — (R' x U)p = 7, consequently, we have 



f{R{x),U{x))dp{x) = / f{R'{x),U{x))dp{x), 
w Jw 

for any bounded, measurable map /oniyx VF. This implies in particular 
f {aoToR) (bo U)dp = ! (a o T o R') (b o [7)cip 

for any bounded measurable functions a and b. Let U' = T o R\ then the above 
expression reads as 

/ a o U b o Udp = / a o U' b o Udp . 
Jw Jw 

Taking a = b, we obtain 

(a o U) (a o [/') dp = ||a o U\\ L 2 {p) \\a o C/'|| L 2 (p) , 



/ 



for any bounded, measurable a. This implies that a o U = a o U' p-almost surely 
for any a, hence U = U' i.e, T o R = To R' p-almost surely. Let us denote 
by S the left inverse of T whose existence follows from Theorem 13.11 and let 
D = {x G W : SoT(x) = x}. Since p(D) = 1 and since R and R' are p-rotations, 
we have also 

P (DnR- l (D)nR'- l (D)) =1. 

Let x e W be any element of D n R~ l (D) n R'~ l (D), then 

= SoToR(x) 
= SoToR'(x) 

= R'(x), 

consequently R = R' on a set of full p-measure.. □ 

Let us give another result of interest as an application of these factorization 
results: it is important to have as much as information about the measures and 
the tranformations which induce them in the setting of Girsanov Theorem, cf. 
and the references there. The problem which we propose is the following: assume 
that, in the case of the Wiener measure, we have a density L with dii(p>, L-p) < 00, 
hence from Theorem 13 .![ a map T : W —>■ W, which is the optimal transport map 
corresponding to the solution of MKP in S(p, L • p). Since the target measure 
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is also spread, the map T possesses a left inverse S such that S o T = Iw p- 
almost surely. Assume now that the transformation T has a Girsanov density, 
i.e., A G L+(/j,), with E[X] = 1 and that 



/ o TA dfx = J f dn , 
for any / G Cb{W). We can now prove: 

Theorem 4.2 Let T be as explained above, assume moreover that 

d H (X ■ /i, /j,) < oo , 

then T has also a right inverse, i.e., T is invertible n- almost everywhere. 

Proof: Denote by G : W — > W the optimal transportation map corresponding to 
the solution of MKP in £(/i, X-fi). Note that both of the measures (T x Iw)(X-fi) 
and (Iw x 0)/x belong to E(/x, A ■ /x). By the uniqueness of the solutions of MKP, 
they are equal, hence, for any a, b G C&(W), we have 



a(T(x))b(x)X(x)dfi(x) = j a(x)b(Q(x))dfi(x) . (4.5) 
Since G(/i) = A • (/,, the equality f)4.5j) can also be written as 

a(ToG(x))b(Q(x))dii(x)= I a(x)b(G(x))dfi(x) . (4.6) 



Since, as T, the map G has also a left inverse, the sigma algebra generated by G 
is equal to the Borel sigma algebra of W, consequently, the relation ()4.6|) implies 
that 

a o T o Q = a , 
/i-almost surely, for any a G C&(W). Therefore we have 

/i ({x G : ToG(x) =x}) = 1, 

since T has already a left inverse, the proof is completed. □ 



4.1 Application to Gaussian measures 

Let us give an example of the above results: Assume that p = /i, i.e., the 
Wiener measure and let K be a Hilbert-Schmidt operator on H . Assume that the 
Carleman-Fredholm determinant det2(/j?+i^) is different than zero, hence the op- 
erator Ih + K : H —>■ H is invertible. Moreover, it follows from the general theory 
that Ih + K has a unique polar decomposition as In + K = (In + K)(Ih + A), 



11 



where In + A is an isometry 3 and In + K is a symmetric, positive operator. 
Note that ^ is compulsorily Hilbert-Schmidt. Let us now define U : W ^ W 
as U(x) = x + <5if(x), where is the if- valued divergence, defined by 

(SK(x),h)H = S(K*h)(x). Then it is known that the measure Ufi is absolutely 
continuous with respect to /z, in fact Ufi is even equivalent to \i since \Ak\ ^ 
//-almost surely, where 

A K = det 2 (I H + K)exp{<5 2 (K) -~\6K\ 2 H j . 

Besides we have 

dU ii 1 
dp, \Ak\ o V ' 

where is the inverse of {7, whose existence follows from the invertibility of 
h->h + 8(K)(x) + Kh on H, cf. |25|. Consequently, 

£[LlogL] = -£[log|A K |]<oc, 

hence dnip, £///) < oo. We shall prove that the polar factorization of U is given 
by 

U=(I W + 5K) o (I w + 5 A) . 
In fact, it follows from Theorem B.6.4 of that 

(I w + 8K) o (I w + 5 A) = I W + 5K + 5A + 5{KA) 

= I W + 5{K + A + KA) 
= I W + 5K. 

Besides V 2 5 2 K = 2K, and since In + K is a positive operator, the Wiener map 
|5 2 /^ is 1-convex, consequently, T = Iw+SK is the transport map and Iw + SA is 
the unique rotation whose existence is proven in Theorem 14.11 The Kantorovitch 
potentials (p and ip of Theorem 13.21 can be chosen as 

<p(x) = ~6 2 K(x) 

for T and 

^x) = - l -5((I H + K)- l K)(x) 

for T" 1 = I w + W- 

Remark 4.1 Let us denote by -Pker<5 the projection operator from ID' (H) to the 
kernel of the divergence operator 5. Then, we have the following identity: 

P kelS (5((I H + K)A)) =6K-6K, 

where K denotes the symmetrization of K. This shows that the polar decompo- 
sition and the Helmholtz decomposition are different in general. 

3 A satisfies the relation A + A* + A* A = 0. 
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We can also calculate the forward Monge-Kantorovitch potential function for the 
singular case as follows: assume that v is a zero mean Gaussian measure on W 
such that <if/(/i, v) < oo. Then there exists a bilinear form q on W* such that 

/ e l ^ x) dv(x) = exp --q(a, a) , 
Jw 2 

for any a G W*. On the other hand, from Theorem 13.21 there exists a (p G D 2 1, 
which is 1-convex, such that Tfi = (Iw + V<p)fi = v. Hence, rewriting the above 
relation with T, we obtain: 

f e i{ta > T{x)) dfx(x) = exp --q(a, a) , (4.7) 
Jw 2 

for any t G H and a G W* . Taking the derivative of both sides twice at t — 0, 
we obtain 

q(a,a) = \a\% + E [(V<p,a)%] +2E[(V<p,ai) H 6a] 

= \a\ 2 H + E [(V(p ®Vtp,a® a) 2 ) + 2E [{V 2 (p, a ® 5) 2 ] , 

where a denotes the image of a under the injection W* •—>■ H. Note that, here, 
V 2 (p is to be interpreted as a distribution. Denote by M the Hilbert-Schmidt 
operator defined by 

M = E [V<p ® Vf] + 2E [VV] • 

We have 

q(a, a) = ((I H + M)a, a) H . 

Let Ih + N be the positive square root of the (positive) operator Ih + M, then 
N is a symmetric Hilbert-Schmidt operator. By the uniqueness of V<£>, it is clear 
that we can choose (p as 

<p = -6 2 N 
* 2 

upto a constant: ip is a 1-convex element of 1)2,1, moreover the map T defined 
by T = Iw + = Jiy + <5iV satisfies the identity ([4.7)1 . hence T is the unique 
solution of the Monge problem and (Iw x T)/x is the unique solution of MKP for 
S(/i,z/). 

5 Strong solutions of the Monge- Ampere equa- 
tion for H- log- concave densities 

Assume that L G ILi x (/x) is of the form 

1 _ f 

L = E\^rf ' 
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where / is an if-convex function. We suppose that L is in ILlogIL(/z), p > 1. 
Denote by <p G D 2 i the potential of the transport problem between fi and v = 
L ■ // which is a 1-convex function. This means that the mapping defined by 
T = Iw + satisfies Tji = L ■ ji and (Iw x T)// is the unique solution of the 
Monge-Kantorovitch problem in z/) with the singular quadratic cost function 
c(x,y) = \x — y\ 2 H . Let A = 1/L o T, it is immediate that A is well-defined since 
LoT ^ /i-almost surely. There exists an inverse map T" 1 , of the form Iw + 77 
defined //-almost everywhere such that T~ x v = fi and that rj G L 2 (u,H). If V 
is closable with respect to u, then g is of the form g = V^> with ip G IL 2 (z/) (cf. 
Remark iH.Hjl . Let L n = £ , [P 1 / n L|'K l ], where 14 is the sigma algebra generated 
by the first n elements of an orthonormal basis (e n ,n > 1) of H and P\/ n is 
the Ornstein-Uhlenbeck semigroup at t — 1/n. It follows from P and from the 
positivity improving property of the Ornstein-Uhlenbeck semigroup, that L n > 
/x-a.s. and that it is of the form \e~^ n i where /„ is a smooth //-convex function 
on W and c = E[e~f}. We denote by (p n , A„, ^ the maps associated to L n , i.e., 
T n = I w + Vy9 n maps \i to the measure L n • \i and 5^ = Iw + V^n maps L n • // 
to /i. Besides, from 0, a; — > x + Vf n (x) is a 1-Lipschitz map, i.e., 

\x + V(f n (x) - y- V(f n (y)\ <\x-y\, 

for any x, y G H™, here it is remarkable that the Lipschitz constant is one and 
it is independent of the dimension of the underlying space. Moreover, we know 
already that the operator valued map Ih + V 2 ya„ > 0, hence <p n is also a concave 
function in D 2i 2- Therefore &p n is a well-defined element of L 2 (/i), |V<p n |g- is 
uniformly exponentially integrable, i.e., there exists some t > such that 

swpE [exp7;|V<^ n |y < 00 , (5.8) 

n 

then the Fatou Lemma implies that 

E [expt\V(p\ 2 H ] 

It follows in particular that (ip n ,n > 1) C ]D Pj2 and it converges to ip in D p l for 
any p > 1, cf., JT]. Moreover, from the Jacobi formula, we have 



/ g o T n A(tp n )dfi = / gd/j, 
Jw Jw 



for any g G Cb(W), where 

A(ip n ) = det 2 ( Ih + VV„) exp |-£(p n - ^\Vf n \ 2 H 

In this case we have A((p n ) = l/L n o T n , and since (L n ,n > 1) is uniformly 
integrable, using the Lusin theorem and the convergence in IL°(//) of (T n ,n > 1) 
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to T, we can show as in Section 5.3 of [26J, that (A(ip n ),n > 1) converges in 
IL O) to 1/LoT. 

The next result about the regularity of <p is of fundamental importance: 

Theorem 5.1 Suppose that L 6 U p> iIL loglL(yu) is H -log- concave. Then the for- 
ward transport potential <p associated to Monge-Kantorovitch problem on L ■ 
fi), belongs to the Sobolev space 102,2- In particular dp £ L 2 {n) and det2(/# + 
V 2 y>) is a well-defined function. 

We shall first prove a lemma of general interest who will imply directly The- 
orem [OJ 

Lemma 5.1 With the hypothesis of Theorem \5.1l it holds that 

£[iV^+||VV|| 2 2] <2£[LlogL], (5.9) 
where \\ ■ \\% denotes the Hilbert- Schmidt norm. 
Proof: Assume first that W = M d and define 
m(t) = -logA(t^) 

t 2 

= tC(p- log det 2 (I + tV i <p) + -\V<p\ 2 H . 
It is immediate to see that 

m'(t) = dp + t trace [(I + tVV)~ 1 (VV) 2 ] + A^^h 

and that 

m"{t) = \\{I + tV\y l {V 2 ip)% + \V V \. 

Note that, since the eigenvalues of V 2 <p are in the interval [—1, 0] , we have m"(t) > 
m"(0). Hence m(l) = m'(0) + |m"(r) > m'(0) + |m"(0), which implies that 

| V<p\ 2 + || VVlla < 2 m(l) = - 2 lo g A (^) - 2 ^ • 
Taking the expectation of both sides gives 

£[|V^| 2 +||VVll2] < -2£[logA^)] (5.10) 

= 2E[L\ogL}. 

□ 

Proof of Theorem I5.lt Recall now that tp n denotes the transport potential 
associated to the measure dv n = L n d[i, where L n = E[Pi/ n L\V n ], hence it follows 
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from the proof of Theorem 4.1 of [llj . that (<p n , n > 1) converges to (p in ID2,i. 
By the Jensen inequality, we get from (|5.10j) 

sup£ [|V<^ n | 2 + HVVJI2] < sup-2E[logA(^ n )] 

n n 

= sup2£[L n logL n ] 

n 

< 2[LlogL], 

consequently p belongs to HD 2 ,2- D 



Corollary 5.1 Under the hypothesis of Theorem \5. 1\ we have 

E[goTA{cp)}<E[g], 

for any positive, measurable function g. In particular A(<p) is a sub-solution of 
the Monge-Ampere equation in the sense that 

A(p) LoT < 1 

pL-almost surely. 

Proof: Since T = Iw + V<p is a monotone shift, the first inequality follows from 
Theorem 6.3.1 of For the second one we use the first one: 

E[goTLoTA(ip)) < E[g L] 

= E[goT], 

for any positive, measurable function g. Since T has a left inverse, the sigma al- 
gebra generated by it is equal to the Borel sigma algebra of W upto /i-negligeable 
sets. Therefore we get the second claim. □ 

We can prove now the main theorem of this section: 

Theorem 5.2 Let LgL logIL(/i) for some p > 1 be given as c _1 e~? , where f is 
an H -convex Wiener function and define the probability measure v as dv = Ldfi, 
where c = E[e~f] is the normalization constant. Let T = Iw + V<£> be the optimal 
transportation of fi to v in the sense of Wasserstein distance, where <p G E^i 
is the forward 1-convex potential function associated to the Monge-Kantorovitch 
problem on v). Then ip G ^2,2 and the Gaussian Jacobian of T exists and 
it satisfies the Monge-Ampere equation: 

k(ip) LoT =1 

ji-almost surely, where 

A(<p) = det 2 (/ H + VV) exp [-Cap - i \Vp\ 2 H \ . (5.11) 
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Proof: We have prepared everything necessary for the proof. First, we can form 
a sequence, denoted by ((p' n , n > 1) such that each tp' n is obtained as a convex 
combination from the elements of the tail sequence (ifk, k > n) and that the 
sequence (p' n ,n > 1) converges to ip in E> 2i2 . Let us denote the Jacobian written 
with ip' n by A n (ip' n ) whose explicit expression is given as 

A{ip> n ) = det 2 (J + VVJ exp |-/y B - ±\Wn\%} 

Let = I w + Vy4 and = I\y + Vt/C Since A — > — logdet 2 (i# + A) is a 
convex function on the space of symmetric Hilbert-Schmidt operators which are 
lower bounded by —In (cf. [3J, p. 63), we have 

- log A(^4) = - log det 2 (l H + *iVVn, 1 

5^ *iV^„ 



< J^-^logA^nJ- 



As we have explained in the paragraph preceding Theorem 15.11 the sequence 
(L n oT n ,n > 1) converges to L o T in probability, and since 

A(ip n ) L n oT n = l 

almost surely, the sequence (— log A((p n ), n > 1) converges to logL o T in prob- 
ability. Moreover, from the construction, (— log A n (ip' n ), n > 1) converges to 
— logA(<£>) in probability. Therefore, it follows from the convexity inequality 
above that 

-logA(p) < logLoT 

almost surely. Consequently (LoT) -1 < A(<£>) almost surely. It follows then from 
Corollary 15. II that (L oT)" 1 = A(<p) almost surely and this completes the proof. 

□ 

The following corollary gives the exact value of the Wasserstein distance: 
Corollary 5.2 With the hypothesis of Theorem \5.<% we have 

\d\{^ L ■ /j,) — E[L log L]+E [logdet 2 (I H + V 2 ^ 
In particular, if 

MA)' 
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where A is an H-convex set, then we have 

Afcp) = ii{A) 

therefore 

1 



fj.(A) = exp t--d A H (ji, L-fi) + £[logdet 2 (/ + V»] 
Proof: Since A = ce^ oT , it follows from the theorem that 

^4(/i,L-/i) = \e[\v<p? h ] 

= -E[foT) -\ogc + E [logdeta^ + VV) 
= E[L log L] + E [log det 2 (/ H + VV)] • 

In particular, the fact that E [logdet 2 (/f/ + V 2 <^)] is always negative explains the 
defect in the Talagrand inequality 23 . To prove the last part, note that we have 
/i(T -1 (v4)) = 1, i.e., 1a ° T = 1 ^-almost surely. Consequently 

1 ■ ^-m 

fi{A) 

Taking the lograrithm of this equality and taking its expectation afterwards gives 
immediately the last formula of the corollary. □ 

Let us give an interesting result about the upper bound of the interpolated 
density whose proof makes use also the convexity results as in the proof of The- 
orem 15.21 : 

Proposition 5.1 Assume the validity of the hypothesis of Theorem \5.S\ suppose 
furthemore that the density L is almost surely bounded, i.e., its exponent is almost 
surely lower bounded by some —a, a > 0. Let T t be defined as T t = I w + tV<p, 
t G [0, 1], then the Radon-Nikodym density L t = d(T t /j,)/d/A is also bounded: 

L t < - exp at 

c 

almost surely, where c = E[exp — /]. 

Proof: Let g be any positive, measurable function on W, by the convexity of 
t — > — log At, we have — log At < — tlogA. Therefore 

E[L t log L t g\ = E[-\ogA t goT t } 

< E[-t\ogAgoT t ] 

= E[-t(foT + \ogc)goT t ] 

< E[(ta- log c)goT t ] 
= E[(ta- log c)L t g] . 
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Consequently 

L t log L t < (ta - log c) L t 
almost surely. □ 



6 Ito-solutions of the Monge- Ampere equation 

In the following calculations we shall take W as the classical Wiener space 
W = C ([0, 1], R), H = H\ i.e., the Sobolev space W 2 ,i([0, 1])- We note that this 
choice does not entail any restriction of generality as indicated in [2H], Chapter 
2.6. Suppose we are given a positive random variables L = whose expecta- 
tion is equal to one, c being the normalization constant. Define the measure v 
as dv = Ldfi. We shall suppose that the Wasserstein distance <i#(/i, v) is finite, 
hence the conclusions of Theorem 13.11 are valid. In order to simplify the discus- 
sion we shall assume that L is strictly positive. The transport map T can be 
represented as T = Iy/ + Vy? again with (p e 1D 21 . Define now 

A = ^. 
LoT 

We have 

J g oT Adfx = J g dfx , 

for any g G (7&(W). This implies that the process (T t ,t e [0,1]) defined on 
[0, 1] x W by 

(t, x) — > T t (x) = x(t) + / D T ip(x)dr , 

Jo 

is a Wiener process under the measure Adfi with respect to its natural filtration 
(J-'f, t G [0,1]), where D t ip represents the Lebesgue density of the map t — > 
W(p(x)(t) G H on [0, 1]. Since T is invertible, we have also 

V ^ = B(W), 

te[o,i] 

upto /x-negligeable sets. Since Adfi is equivalent to the Wiener measure, the 
process (T t ,t G [0, 1]) is a /z-semimartingale with respect to its natural filtration. 
It is clear that it has a decomposition of the form 

T t = Bj + A t , 

with respect to /i, where B T is a /i-Brownian motion and A is a process of finite 
variation. Since we are dealing with the Brownian nitrations, (A t ,t G [0,1]) 
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should be absolutely continuous with respect to the Lebesgue measure dt of [0, 1]. 
In order to calculate its density it suffices to calculate the limit 



lu.- h E[T t+h -T t \?n- 
To calculate this limit, it is enough to test it on the functions of the type g o T t : 
E[(T t+h -T t )goT t ] 



E[(W t+h -W t )goW t L] 
E[(6U [t:t+h] )goW t L] 
E[(U [t , t+h] ,V(LgoW t )) H ] 



E 



goW t 



t+h 



D T Ldr 



(6.12) 
(6.13) 



where U[ t> t+h] is the element of H whose Lebesgue density is equal to the indicator 
function of the interval [t, t + h]. Note that for the equality (|6.12|) . we have used 
the fact that 5 = V* under the Wiener measure fi and the equality ()6.13|) follows 
from the fact that the support of \7(g(W t )) lies in the interval [0,t], hence its 
scalar product in H with U{ t ,t+h\ is zero (cf.|25|). Hence we have 

\iw^E[T t+h -T t \Fj] = -E[D t foT\rJ] 



-E v [D t f\T t ]oT, 



dt xd/x-almost surely, where the last inequality follows from the fact that T 1 (jF t ) = 
Tj. Hence we have proven 

Proposition 6.1 The transport process (T t ,t e [0, 1]) is a (/a, {Tf))-semimartingale 
with its canonical decomposition 

T t = Bf-f E v [D T f\T T ]oT dr 
Jo 

= Bj- [ E[D T foT\F?] dr. 
Jo 

We can give now the Ito solution of the Monge- Ampere equation: 

Theorem 6.1 Assume that f e 1)2,1 be such that c = E[exp(—f)] < oo, denote 
by L the probability density defined by -e~~' and by v the probability dv = Ldji. 
Assume that c£ff(//, v) < 00 and let T = I\y + Vy? be the transport map whose 
properties are announced in Theorem \3. 6 A We have then 

A = exp|^ E u [D t f\T t ]oTdBj - - jf E u [D t f\F t ] 2 oT dt\ . (6.14) 
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Proof: From the Ito representation formula [21] > we have 

L = exp j-^ E v \D t f\T t \dW t - 1 -^ E v {D t f\T t f dt^ . 
Since the Girsanov measure for T has the density A given by 

a = -!-, 

LoT 

we have, using the identity T~ l {!Ft) = and Proposition EH1 

LoT = expj- J^E v [D t f\T t ]oTdT t -^j^E v [D t f\T t \ 2 oTdt^ 
= ex p{-_^ E u [D t f\T t ]oT (dBj - E u [D t f\T t }oTdt) 
jf E v [D t f\F t foTdt 



which is exactly the inverse of the expression given by the relation (|bM4j) . □ 

The following proposition explains the relation between the semimartingale rep- 
resentation of T and the polar factorization studied in Section 

Proposition 6.2 Let X be the process defined by 

X t = W t + [ E u [D T f\F T }dr, 
Jo 

then ToX is a v -rotation, i.e., ToX(u) = v, in fact it is the minimal v -rotation 
in the sense that 

inf E v [\0 - X\l] = E U [\T o X - X\l] , 

where 1Z V denotes the set of transformations preserving the measure v . Finally 
the Brownian motion B T is the rotation corresponding to X o T. 

Proof: Since E[L] = 1, v is the Girsanov measure for the transformation X, 
consequently, we have 

E v [g(ToX)} = E[g(ToX)L] 

= E[g(T)) 
= E u [g] } 
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for any g G CbiW) and this implies T o X(u) = v. Let now O G 1Z U , then the 
measure O x X(u) belongs to fi). Since T x !(//) is the solution of MKP in 
£(z/, /i), we have 

E v [\0 - X\%] > E V [\T o X - = d H (ji, vf . 

The uniqueness follows from the same argument as used in the proof of Theorem 
14.11 The last claim is obvious since X o T is a //-rotation, hence as a process it is 
a Brownian motion, then by comparing it with the result of Proposition 16. 1[ we 
see that B T = X o T. □ 

Let us give some immediate consequences of these results whose proofs follow 
from the results of this section and from Theorem 15.21 : 

Corollary 6.1 We have the following identity 



■\ogE[e~ f ] = E 
= E 



foT + \ [ E u [D t f\F t } 2 oTdt 
1 Jo 

foT + ±^E[D t foT\fT] 2 dt 



If, furthermore, f is H-convex, then we also have 



-log£[e -/ ] = E 



1, 



foT- \ogdet 2 (I H + VV) + j I V ^Ih 



In particular we have the exact characterization of the Wasserstein distance be- 
tween /j, and v: 



X -d\{^v) = £[logdet 2 (itf + VV)]+^ 



E v \D t f\T t \ 2 oTdt 



E [logdet 2 (/H + VV)] + £[£logL] 
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